Topological Crystalline Insulators in the SnTe Material Class 



Timothy H. Hsieh, 1 Hsin Lin, 2 Junwei Liu, 1,3 Wenhui Duan, 3 Arun Bansil, 2 and Liang Fu 1 '^ 

department of Physics, Massachusetts Institute of Technology, Cambridge, MA 02139 
2 Department of Physics, Northeastern University, Boston, Massachusetts 02115, USA 
3 Department of Physics and State Key Laboratory of Low-Dimensional 
Quantum Physics, Tsinghua University, Beijing, 10084, China 

Topological crystalline insulators are new states of matter in which the topological nature of 
electronic structures arises from crystal symmetries. We theoretically demonstrate that SnTe is a 
topological crystalline insulator due to its intrinsically inverted band structure at four L points in 
the Brillouin zone. As a consequence, we predict that SnTe has metallic surface states with an even 
number of Dirac cones on crystal surfaces such as [001] or [111], which are symmetric about the 
[110] mirror plane. These surface states, exhibited by our first-principles calculations, form a new 
type of Dirac metal protected by mirror symmetry and robust against disorder. We also find that 
the related semiconductors PbTe and PbSe become topological crystalline insulators under pressure 
after their band gaps are inverted. 

PACS numbers: 



The discovery of topological insulators [THS] has at- 
tracted much interest in the topological classification of 
electronic band structures beyond the existing Z 2 mate- 
rial class. In recent years, theoretical explorations con- 
tinue to find new phases of matter protected by symme- 
try operations other than time reversal transformation. 
A fruitful direction is to incorporate crystal symmetries 
into the classification scheme. It has been recognized 
that crystal symmetries play a multifaceted role: they 
can constrain the band topology@HS] or engender new 
classifications [7J [8]. The latter gives birth to new topolog- 
ical phases dubbed topological crystalline insulators [5], 
which have metallic surface states protected by crystal 
symmetries. Thanks to the richness of crystal structures, 
a large number of topological crystalline insulators awaits 
discovery — both theoretically and experimentally. 

In this Letter, we theoretically demonstrate that SnTe 
is a topological crystalline insulator with mirror symme- 
try. As a consequence, it is predicted to have robust sur- 
face states with an even number of Dirac cones on crystal 
surfaces such as [001] or [111], which are symmetric about 
any of the six equivalent [110] mirror planes. The exis- 
tence of these surface states is dictated by a nonzero in- 
teger topological invariant — the mirror Chern number [S]. 

SnTe has a simple rocksalt structure (Fig. la); its fun- 
damental band gaps are located at four equivalent L 
points in the face-centered-cubic (FCC) Brillouin zone. 
It has long been established that the ordering of the con- 
duction and valence bands at L points in SnTe is in- 
verted relative to PbTe, so that the band gap of the alloy 
Pbi-^Sn^Te closes and re-opens as x increases [TO]. Since 
this band inversion occurs at an even number of points, 
neither SnTe nor PbTe in the rocksalt structure is a topo- 
logical insulator with Z 2 topological order. 

However, the above band inversion has significant con- 
sequences when the mirror symmetry of the FCC lattice 
is taken into consideration. Consider the plane YL\L 2 in 




FIG. 1: (a) the crystal structure of SnTe; (b) the FCC Bril- 
louin zone showing the plane TL1L2 which is invariant under 
reflection about the (110) axis and projects onto the TXi line 
in the [001] surface. 



momentum space, defined by the three points r,Li and 
L 2 (Fig. lb). Crystal momenta on this plane are invari- 
ant under reflection about the [110] mirror plane in real 
space. This allows us to label the Bloch wavefunctions on 
this plane by their eigenvalues ±i under the mirror opera- 
tion M (which satisfies M 2 = — 1 for spin 1/2 electrons). 
Each class of M = ±i mirror eigenstates has an associ- 
ated Chern number n±i, and the mirror Chern number 
^m[S] is defined by tin = {n+i — n^i)/2. Provided that 
mirror symmetry is present, um is an integer topologi- 
cal invariant. A nonzero mirror Chern number defines a 
topological crystalline insulator with mirror symmetry. 

We first demonstrate that SnTe and PbTe have mirror 
Chern numbers which differ by two, and hence, one of 
them is a topological crystalline insulator. This is estab- 
lished by considering the band inversion between SnTe 
and PbTe, which is captured by the k ■ p theory at L 
points [TT]: 

H = ma z + v(k 1 s 2 - k 2 s 1 )a x + v 3 k 3 a y . (1) 
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Here k\ , k-2 , k 3 form an orthogonal system with k 3 along 
TL and fci along (110); a z = ±1 corresponds to the 
p-orbital on the cation (Sn or Pb) and anion (Te) re- 
spectively; s 3 — ±1 labels the total angular momentum 
j = ±1/2 along TL. A positive to means that the con- 
duction and valence bands at L are respectively derived 
from the cation and anion, and vice versa for negative 
to. The form of H in Eq.fJTJ) is uniquely determined 
by the point group symmetries which leave an L 
point invariant [1 2 . In particular, reflection about the 
[110] mirror plane is represented by M = —is\. On the 
mirror-invariant plane TL1L2 (k\ — 0), H reduces to 
Hq = m<7 z — vk 2 s\o x + v 3 k 3 a y . Due to mirror symmetry, 
Hq decomposes into the s% = 1 (M = —i) and Si = — 1 
(M = i) subspaces: 

H = ma z =F vk 2 a x + v 3 k 3 a y , (2) 

each of which describes a 2D massive Dirac fermion. 

In going from PbTc to SnTe (increasing x in 
Pb 1 _ ;r Sn :E Te), the cation/anion character of the conduc- 
tion/valence bands becomes switched at L, which in the 
k ■ p theory corresponds tom-> — m. The sign reversal 
of to at one L point changes the Chern number of the 
s% = ±1 subspace by ±1, and hence the mirror Chern 
number changes by one. Furthermore, L\ and L 2 are 
related by a two-fold rotation around the (110) axis. Un- 
der this rotation, both the spin polarization s\ and the 
orientation of the TLiL 2 plane are flipped. As a result, 
the Berry curvatures at L\ and L 2 are related by 

F +i {L x ) = -F-i{L 2 ), F-iiLx) = -T +l (L 2 ). (3) 

so that 

( T +i - 7-_0k = (J r +i -^-i)|i a , (4) 

i.e., both Li and L 2 contribute equally to the change in 
mirror Chern number. The net result is that the band in- 
version changes the mirror Chern number for the TL\L 2 
plane by two. This means that in the presence of the 
[110] mirror symmetry, SnTe and PbTe are topologically 
distinct. 

In order to determine which of the two is topo- 
logically nontrivial, we analyze their band structures 
obtained from first-principles density functional the- 
ory calculations [T3HH] using the generalized gradient 
approximation [IB]. The results we obtain are consistent 
with several previous studies [T?H2T)] . In particular, it is 
generally agreed that the band ordering at L point is 
correct at the total-energy optimized volume. 

Based on two first-principles findings, we now demon- 
strate that, of the two materials, SnTe is the topologically 
nontrivial one. First, the conduction and valence bands 
of PbTe throughout the Brillouin zone are primarily de- 
rived from the p-orbitals of Pb and Te atoms respectively 
(Fig. 2b). This suggests that PbTe is smoothly connected 
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FIG. 2: Electronic structures and band-gap evolution of 
SnTe and PbTe: The electronic band structures of SnTe and 
PbTe are shown in a and b, respectively. The size of the red 
dots represents the fraction of electronic charge residing on 
Te atoms. The exchange of the band character at L point 
as highlighted in the grey area indicates the intrinsic band 
inversion of SnTe. The band gap energy as a function of 
the lattice constants is shown in c. The negative gap area 
indicates the topological crystalline insulator phase. 



to the atomic limit, in which Pb orbitals are empty and 
Te orbitals filled due to their on-site energy difference. 
In contrast, the orbitals in the band structure of SnTe 
are switched near L points: the conduction band edge 
is derived from Te and the valence band edge from Sn 
(Fig. 2a). Therefore SnTe has an intrinsically inverted 
band structure (to < 0). Our conclusion is further sup- 
ported by the dependence of the band gap on the lattice 
constant obtained in first-principles calculations (Fig. 2c) 
and measured via pressure coefficients [UJ. As the lat- 
tice constant increases, the band gap of PbTe increases 
monotonically, whereas that of SnTe decreases to zero 
and then re-opens. This gap closing signals a topological 
phase transition, in which SnTe at ambient pressure is a 
topological crystalline insulator with the aforementioned 
mirror Chern number um = 2. 

The nonzero n^j in SnTe dictates the existence of sur- 
face states on any crystal surface which is symmetric 
about the [110] mirror plane (or its equivalents). Two 
common surface terminations satisfying this condition 
are [001] and [111]. Using bulk-boundary correspon- 
dence, we now infer the topology of these surface bands. 
For the [001] surface, the plane YL\L 2 in the bulk Bril- 
louin zone projects onto the symmetry line YX\ in the 
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surface Brillouin zone, with both L\ and L2 projecting 
onto Xi . riM = 2 dictates that along the line X\ — T — X\ 
there must exist two pairs of counter-propagating chiral 
surface states with opposite mirror eigenvalues. By rota- 
tional symmetry, such surface modes also appear along 
the line X2 — f — X2. But they are absent along any 
other mirror-invariant line. Since the crossing of two 
mirror branches forms a Dirac point, the [001] surface 
states have four Dirac points on the four equivalent FX 
lines. Similar considerations apply to the [111] surface. 
For this surface, the TLiL 2 plane projects onto the high- 
symmetry line TM. Therefore two surface modes exist 
along the three equivalent lines M — T — M. In contrast 
to Z2 topological insulators, the Dirac points here are lo- 
cated at an even number of non-time-reversal-invariant 
momenta. Our theoretical prediction of these surface 
states in SnTe is the main result of this work. 

The above surface states are explicitly demonstrated in 
a slab geometry along the (001) axis using first-principles 
calculations. As predicted by the mirror Chern number 
jim = 2, two surface bands with opposite mirror eigen- 
values cross each other to form a Dirac point on the line 
r and X (Fig. 3a). The Dirac velocity is found to be 
1.7 x 10 5 m/s in the TX direction. Interestingly, these 
surface states exhibit a Lifshitz transition — a change of 
Fermi surface topology as a function of Fermi energy. 
Fig. 3b shows a set of Fermi surfaces at different energies. 
As the Fermi energy is lowered from the Dirac point, the 
Fermi surface near X undergoes a Lifshitz transition: two 
disconnected Fermi surfaces outside X touch each other 
and reconnect to form one large and one small pocket, 
both enclosing X. 

In order to understand both the connection between 
the bulk and surface bands and the effect of potential 
perturbations at a microscopic level, we introduce a sim- 
plified tight-binding model for SnTe. This model is con- 
structed from the Wannier functions of the conduction 
and valence bands, which are primarily three p-orbitals 
of Sn and Te atoms. The Hamiltonian H t b is given by 

j r,a 

(r') + h.c. 

j,j' (r,r').a 
j r,a,/3 

Here r labels the site, j = 1,2 labels the Sn or Te atom, 
a =t, l labels electron's spin. d rr ' is the unit vector con- 
necting site r to r'. The components of vectors and 
c correspond to the three p-orbitals. In the Hamiltonian 
([5]), m is the on-site potential difference between Sn and 
Te; ti2 — <2i is the nearest-neighbor hopping amplitude 
between Sn and Te; tn and £22 are the next nearest- 
neighbor hopping amplitudes within a sublattice; the Ai_2 
term is L-s atomic spin-orbit coupling. The bulk and sur- 
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FIG. 3: The [001] surface states of SnTe: (a) band disper- 
sion and (b) Fermi surface. Note that in between F and X 
two surface bands with opposite mirror eigenvalues cross the 
Fermi energy, in agreement with um = 2. (c) a set of Fermi 
surfaces at different energies, exhibiting a Lifshitz transition. 

face bands of the above tight-binding Hamiltonian nicely 
reproduce the essential features of the first-principles cal- 
culation. 

Using this tight-binding model, we now study the elec- 
tronic properties of the (001) surface states under var- 
ious perturbations. The four Dirac points located at 
k J , j = 1, ...4 are related by the four-fold rotation C4. 
For convenience, we choose a basis for each Dirac dou- 
blet ip a (k 3 ), a = ± which satisfies 

C 4 : MW) -> e- M7r/4 V Q (k J+1 ). (6) 

In this basis, the k ■ p Hamiltonians at four Dirac points 
take an identical form with the same chirality: 

H S { = v±kis 2 - v\\k 2 si, (7) 

Here k\ and k 2 form a local right-handed coordinate sys- 
tem centered at each k J , with fc 2 parallel to k 7 . Provided 
mirror symmetry is present (fci — > — fci,s 2 — > — S2), the 
Dirac points can move along the TX line, but cannot 
annihilate with each other. 

Perturbations which break mirror symmetries can gen- 
erate mass terms rrijS^ at the four Dirac points k 7 and 
open up gaps. The nature of the resulting gapped phase 
is determined by the relative signs of rrij, which in turn 
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FIG. 4: Under a ferroelectric distortion with a displacement 
u, the four surface Dirac points (detailed on the right) are 
gapped with mass signs depending on the direction of uy. 
This leads to the four gapped phases (I-IV). 

depend on the symmetry of the perturbation. For exam- 
ple, a perpendicular Zccman field B is rotationally in- 
variant, and hence yields m\ — mi = 777,3 — m 4 B. 
This drives the surface into an integer quantum Hall 
state. Another example is a structure distortion in which 
atoms are displaced by u. In fact, SnTe distorts along 
the (111) direction into a rhombohedral structure at low 
temperature. This perturbation is time-reversal invariant 
but breaks rotational symmetry, leading to a ferroelectric 
phase. The effect of the distortion is captured by adding 
to the tight-binding model ([5| a modulated hopping term 

H^H + 5t c L( r )"drr' u-c 2a (r') + h.c, (8) 

(r,r'),a 

By symmetry analysis ([6| , the 5t term gives the following 
mass terms to the surface states: 

777 j oc (u x kg) • z, (9) 

where z is the surface normal. It follows that for the 
rhombohedral distortion u cx (111), the two Dirac points 
along (110) direction are gapped, but the other two along 
(110) remain gapless[2"2"]. 

In general, all four Dirac points arc gapped by the dis- 
tortion unless its in-plane component un coincides with 
the direction of any kg. Therefore, there are four types 
of fully-gapped surface ferroelectric phases depending on 
the direction of um (Fig. 4). The domain wall between 
two such phases (A and B) is particularly interesting, if 
ujj 1 belong to the I-III quadrants and Uj^ the II-IV quad- 
rants. This domain wall carries two branches of counter- 
propagating spin-polarized electrons: each branch arises 
from the sign reversal of Dirac mass across the domain 
wall[23]- Such one-dimensional electrons are protected 
against backscattering by time reversal symmetry |24j. 
and hence form perfectly conducting channels. 

The above analysis of surface states under mirror sym- 
metry breaking also provides insight into their robustness 
against time-reversal-invariant disorder. By definition, 
disorder breaks mirror symmetry in a random way lo- 



cally, not macroscopically. Let us assume that the dis- 
order potential is slowly varying. The surface is then 
an equal-weight mixture of the type I-III and type II- 
IV ferroelectric domains mentioned above, with domain 
walls percolating throughout the surface. Since the do- 
main wall is a perfectly conducting channel, we expect 
that the entire surface is conducting. The above rea- 
soning resembles recent studies of disordered weak topo- 
logical insulators|i25-28 . Furthermore, quantum tunnel- 
ing between the domain walls can be accounted for in a 
Chalker-Coddington type network model, which will be 
studied in the future. 

On the experimental side, the surface states of SnTe we 
predicted can be readily detected in angle-resolved pho- 
toemission spectroscopy (ARPES) and tunneling spec- 
troscopy experiments. We note that a photoemission 
experiment 29 was recently performed on SnTe. The un- 
derlying mirror Chern number can be deduced from the 
spin polarization of surface states as in Sb[30]. More- 
over, SnTe-based thin films and superlattices have re- 
markably high mobility exceeding 2500cm 2 /Vs at room 
temperature|31j . which provide a promising platform for 
device applications. 

Finally, we relate our work to a wider class of materials, 
including PbTe and PbSe. While both are topologically 
trivial at ambient pressure, our first-principles calcula- 
tion (Fig. 2c) shows that decreasing the lattice constant 
by 2% inverts the band gap and drives PbTe into a topo- 
logical crystalline insulator. This band inversion is real- 
ized under moderate pressure (around 3GPa in PbTe and 
2GPa in PbSe OET]). Alternatively, one can achieve 
the topological regime by growing these materials on sub- 
strates with smaller lattice constants. As a precedent, 
high-quality PbTe quantum wells have been fabricated 
and exhibit ballistic transport [321 133) . 

We briefly comment on how our work relates to early 
pioneering works, which predicted the existence of two 
dimensional massless Dirac fermions at the interface of 
PbTe and SnTe[Ml 133, or domain wall of PbTe[H]. 
Our work has made it clear that only interfaces sym- 
metric about the (110) mirror plane have protected gap- 
less states, which are solely derived from the topological 
crystalline insulator SnTe and exist even when PbTe is 
removed. In light of their topological nature which we 
identified, SnTe and pressurized PbTe are likely to lead 
a new generation of topological materials. 
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